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Abstract. A Key feature of the minimal supersymmetric extension of the Standard Model 
(mssm) is the existence of a light Higgs boson, the mass of which is not a free parameter but 
an observable that can be predicted from the theory. Given that the lhc is able to measure 
the mass of a light Higgs with very good accuracy, a lot of effort has been put into a precise 
theoretical prediction. 

We present a calculation of the SUSY-QCD corrections to this observable to three-loop order. 
We perform multiple asymptotic expansions in order to deal with the multi-scale three-loop 
diagrams, making heavy use of computer algebra and keeping a keen eye on the numerical error 
introduced. 

We provide a computer code in the form of a Mathematica package that combines our three- 
loop SUSY-QCD calculation with the literature of one- and two-loop corrections to the Higgs mass, 
providing a state-of-the-art prediction for this important observable. 



1. Introduction 

The minimal supersymmetric extension of the Standard Model (mssm) is a promising candidate 
for physics beyond the Standard Model. Its Higgs sector is a two-Higgs doublet model with 
the additional constraint that supersymmetry relates the quartic Higgs couplings to the gauge 
couplings of the theory. This increases the predictiveness of the model and allows the Higgs 
sector to be parametrised by just two new parameters, the mass Ma of the pseudoscalar Higgs 
and the ratio tan f3 of the vacuum expectation values of the Higgs doublets. 

In particular, the mass Mh of the light Higgs boson is not a free parameter, but can be 
predicted. At tree level, Mh is bounded above by Mz, but radiative corrections shift the value 
significantly. Since Mh will be a precision observable once the Higgs is found at the Large 
Hadron Collider (lhc), it is imperative to have a precise theoretical prediction. Consequently, 
a lot of effort has been put into the calculation of radiative corrections to Mh at the one- and 
two- loop level (see, for example [HEllSllllElEllTllH]). The remaining theoretical uncertainty 
has been estimated to be about 3 — 5 GeV [21 [TU] , which is confirmed by a study of the leading 
and next-to-leading terms in ln{MsusY /Mt), where MsusY is the typical scale of SUSY particle 
masses, at three- loop order 

This uncertainty of the theoretical prediction justifies a calculation of the next term in the 
perturbative expansion. A study of the corrections to Mh shows that the contributions from 
loops of top quarks and their superpartners, the stops, are dominant at the one- and two-loop 
level. In |l2l[T3], we have calculated three-loop SUSY-QCD corrections to these diagrams. The 
calculation of these terms is the subject of this talk. 



2. Calculation of the Three-Loop Corrections 

Motivated by the discussion above, we restrict the calculation at the three-loop level to the 
SUSY-QCD corrections where the Higgs couples to a top quark or its superpartners. In the 
perturbative expansion of M^, these are the terms of O [ata^), where Ug is the strong coupling 
and at denotes the coupling of the Higgs to the top quark. 

We thus have to evaluate three-loop propagator diagrams in SUSY-QCD, which faces us with 
two complications. First, we need a convenient regulator that respects supersymmetry, and 
second, a lot of masses appear in our diagrams. 

2.1. Regularisation by Dimensional Reduction 

The regularisation workhorse of multi-loop calculations is dimensional regularisation (dreg) [T^, 
where the number of spacetime dimensions is altered from 4 to d = 4 — 2e, and the divergences 
of the loop integrals are manifest as poles in e. 

Unfortunately, DREG does not respect supersymmetry. An easy way to see this is that 
changing the number of spacetime dimensions also changes the number of degrees of freedoms of 
the vector fields, and supersymmetry requires an equal number of bosonic and fermionic degrees 
of freedom. While it is possible (but tedious) to manually construct finite counterterms that 
restore the supersymmetric ward identities, a more convenient approach was suggested by Siegel 
under the name of dimensional reduction (dred) |15j . The main idea is that the change of 
dimensions need only affect momenta in order to regularise the loop integrals. In dred, one 
splits the four dimensional space into orthogonal spaces of dimension 2e and d = 4 — 2e, and 
restricts the loop momenta to the d dimensional subspace while all gauge fields are kept four 
dimensional. The 2e components of the gauge fields transform as a 2e tuple of scalar fields and 
are called e-scalars. 

The consistent formulation of DRED is more involved than sketched here [T6] , and the question 
whether the inclusion of the e-scalars truly restores the supersymmetric ward identities in all 
cases is not yet resolved. However, DRED has successfully been applied in many multi-loop 
calculations [13 [El [201 [IB [22l [231 [211 [251 [26] , so far without observing a violation of the 
ward identities. 

2.2. Asymptotic Expansion in the Masses 

The difficulty of calculating a Feynman diagram rises with the number of loops, the number of 
external momenta and the number of masses in the diagram. To obtain the O [atOd'^) terms of 
Mfi, we have to evaluate diagrams with three loops, two legs and a lot of masses: treating the 
light quarks as massless and their superpartners as mass degenerate with mass Mq, we are left 
with five masses: the top quark mass Mt, the masses of its superpartners M^- , Mf^, the mass of 
the gluino Mg and Mg. Calculating these diagrams without any approximations is not feasible 
with contemporary methods. 

Given the mild dependence of the one-loop corrections to M/j on the external momentum, 
it is reasonable to expand the diagrams in small external momentum, reducing the problem to 
the evaluation of vacuum diagrams. A further simplification is possible by performing nested 
asymptotic expansions [271 |28l [29l |30j , which expresses the multi-mass diagrams as a series in 
ratios of the masses and logarithms of these ratios. The coefficients in the series involve only 
one-scale integrals. Assuming that the ratios of the masses are small, the series should give good 
approximations to the original diagrams. 

Of course, before carrying out the asymptotic expansions, one has to identify which mass 

ratios are small. Given that so far none of the superpartner masses are known, this is an 

undecidable problem. So, instead of calculating with a fixed, known, hierarchy among the 

masses, one has to systematically consider various possible mass hierarchies and carry out the 

calculation for each of these. To get a prediction of Mh for a specific point in the parameter 
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Figure 1. Behaviour of the perturbative series for Mh using the on-shell and the minimal 
subtraction scheme (for degenerate superpartner masses). The short-dashed, long-dashed and 
solid lines are the one-, two-, and three-loop approximation to Mh, respectively. In the on-shell 
scheme, the distance between one- and two-loop approximation is much larger than in the br 
scheme. 



space of the mssm, where the masses of the superpartners have specific values, one has to choose 
whichever hierarchy fits these values best. 

We carry out the calculation using the following setup: In a first step, all diagrams 
contributing to at O (ata^) are found using the program QGRAF [31]. There are 30.717 
diagrams. For each hierarchy, these diagrams are expanded asymptotically using q2e/exp [32]. 
The one-scale integrals are then calculated using the form |33j program matad [34j . 



2.3. Renormalisation 

For the renormalisation of the parameters entering our calculation, we adopt a variation of the 
i5R-scheme, i.e. minimal subtraction using dimensional reduction. This leads to a much better 
behaviour of the perturbative series compared to using on-shell renormalisation (see Fig. [T]). 



3. Computer Code 

In order to provide a precise prediction for the value of within the mssm, our O (ata^) 
terms have to be combined with the wealth of corrections from other sectors of the MSSM 
at one- and two-loop level that are available in the literature [HEIEIIIIEIEIITIIH]. Also, 
the choice of hierarchy mentioned in 2.2 should be automatised. In p!3], we have presented 
the Mathematica package H3m, which is publicly available, to address these points. For 
convenience, it implements the Susy Les Houches Accord SLHA [35] and has an interface to the 
spectrum generators SOFTSUSY [36], SuSpect [37] and SPheno [38]. 

The program uses FeynHiggs |39l [9] to get a two- loop approximation of Mh. An obstacle 
for adding our terms of O (^atOg) to the result delivered by FeynHiggs is the usage of on-shell 
renormalisation in FeynHiggs. In order to be consistent, we have to convert the O (at) and 
O (atas) terms that contain on-shell parameters to the dr scheme. Thanks to [7], which gives 
a compact expression for these terms both in the dr and in the on-shell scheme, this obstacle is 
easily overcome. 

The corrections to Mh depend very strongly on the mass Mt of the top quark and the strong 
coupling as- These parameters must be known, within SUSY-qcd renormalised in the i5r scheme, 
as precisely as possible. In |30], the relation between the top mass in the i5r- and the on-shell 
scheme in SUSY-qcd has been calculated to O (a^). Using the library tsil [41J, this relation 

can be used to obtain MJ^^ from the known value of the top quark pole mass. 

We determine ag in SUSY-qcd from the value of as in five-flavour QCD at the Z mass 
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following |42| by running, within five-flavour QCD, to the decoupling scale where we perform 
the transition to the full theory. We then run, within SUSY-qcd, to the desired value of the 
renormalisation scale. 

In order to choose the most appropriate hierarchy and get an estimate for the error introduced 
by the asymptotic expansion, we compute an approximation to the two-loop corrections of 
O (atas) using asymptotic expansions in the mass ratios and compare this approximation to 
the result of [7|. The hierarchy that matches the exact two-loop result best is chosen for the 
calculation of the three-loop term, and the error at two-loop level is recorded to get a handle on 
the error in the three- loop result (see Fig. |4|. 

4. Numerical Results 

In this chapter, we present a state-of-the-art numerical prediction for Mh in the MSUGRA model 
that has been obtained with H3m. Fig. l2|shows the dependence of Mh on the parameters mi/2, 
mo, tan/3, and for fi> 0. The shadea bands around the individual curve show the variation 
of Mh when Alt is varied between 171.4 GeV and 174.4 GeV. We restrict the plot to values of 
mi/2 > 300 GeV in light of the latest exclusions from the lhc experiments [ISl HH Ii5l H6] . 

We observe that of all MSUGRA parameters, varying mi/2 has by far the largest impact on 
the Higgs mass. We also observe that the present uncertainty of the top mass directly translates 
to a parametric uncertainty of about one GeV of the Higgs mass. 

To estimate the effect of unknown higher orders. Fig. [3] shows the two- and three-loop 
corrections to Mh- We observe a partial cancellation between the two- and three- loop correction. 
We also observe that while the higher order corrections do decrease in magnitude, they do not 
do so by a large factor. Thus, we should be careful when estimating the magnitude of the 
missing higher order contributions. We choose to assign 50 % of the three-loop contribution as 
a theoretical error. This amounts to a hundred MeV for low values of mi/2 and to about one 
GeV for mi/2 above one GeV. Using the scale variation as an error estimate would lead to a 
smaller error |13| . 

By performing asymptotic expansions, we have introduced an additional source of uncertainty 
to the prediction of Mh- With Fig. [4j we can analyse how large this uncertainty is. It shows 
the error that we would have introduced had we made the same approximation at the two-loop 
level. We see that this error is typically at or below 100 MeV, or 200 MeV for low values of mi/2 
and tan/3. Since the three-loop corrections are smaller than the two-loop corrections, we expect 
the actual error due to asymptotic expansions to be below 100 MeV. 

5. Conclusions 

We have presented a calculation of the O [ata'^) corrections to the mass Mh of the light Higgs 
boson in the mssm. These contributions shift the value of Mh by 1 — 3 GeV, depending on the 
mass spectrum of the superpartners. 

The results have been implemented in the program H3m, which is freely available [13]. 

Using our results, we improve the theoretical error significantly. The theoretical error is now 
on the order of about 100 MeV for light and 1 GeV for heavy superpartner masses. This is 
comparable to the parametric uncertainty with the top mass and a^. 
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Figure 2. Theoretical prediction for the light Higgs mass in the MSUGRA model. We plot Mh 
over the fermion mass parameter mi/2- The panels show, from left to right, an increasing value 
of tan/3, and, from top to bottom, of the scalar mass parameter mo- We show three curves 
according to different values of the trilinear coupling Aq. 
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Figure 3. Two- loop (dashed) and three- loop (solid) contributions to M^. The same conventions 
as in Fig.[2]are adopted. For lucidity, only the curves for Aq = and Mf = 172.9 GeV are shown. 
We observe that the three-loop terms are smaller than the two-loop terms, but not by a large 
factor. This suggests a conservative approach when estimating higher order effects. 




Figure 4. Deviation of the expansion in masses from the exact result at two loops. The 
conventions are the same as in Fig. [2j The error introduced by the asymptotic expansions is 
well under control. The little bumps in the curves are due to a change from one hierarchy to 
another. 
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